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SUMMERY 

In  this  paper,  we  consider  a  low-complexity  method 
of  weighted  subspace  fitting  (WSF)  for  direction-of-arrival 
(DOA)  estimation.  With  the  properties  of  the  multi-stage 
wiener  filter  (MSWF),  we  derive  a  novel  criterion  function 
for  the  WSF  method  without  the  estimate  of  an  array  covari¬ 
ance  matrix  and  its  eigendecomposition.  A  new  approach  for 
noise  variance  estimation  is  also  proposed.  Numerical  results 
indicate  that  by  selecting  a  specific  weighting  matrix,  the 
low-complexity  WSF  estimator  can  provide  the  comparable 
estimation  performance  with  the  conventional  WSF  method. 

I.  INTRODUCTION 

Super-resolving  correlated  or  even  coherent  signals  is 
the  fundamental  problem  in  array  signal  processing,  which  is 
frequently  encountered  in  many  areas  such  as  communication, 
radar,  sonar  and  geophysical  seismology.  The  subspace  based 
methods,  which  resort  to  the  decomposition  of  the  observation 
space  into  signal  subspace  and  noise  subspace,  can  provide 
the  outstanding  estimation  performance.  In  the  literature,  the 
classical  subspace  based  methods  have  been  investigated  ex¬ 
tensively  [  1  ]-[2] .  Nevertheless,  the  classical  subspace  based 
methods  involve  the  estimate  of  an  array  covariance  matrix  and 
its  eigenvalue  decomposition  (EVD),  which  is  rather  computa¬ 
tionally  intensive  for  the  case  where  the  model  orders  in  these 
matrices  are  large.  To  reduce  the  computational  complexity 
of  the  classical  subspace  based  methods,  a  number  of  low- 
complexity  methods  without  eigendecomposition  have  been 
proposed  in  [3]-[4],  Normally,  the  existing  linear  operation 
based  methods  with  low  complexity  [3] -[6]  find  the  signal 
or  noise  subspace  by  a  partition  of  array  response  matrix 
or  exploiting  the  array  geometry  and  its  shift  invariance 
property  [6],  and  then  estimate  the  directions  of  arrival  (DOAs) 
of  signals  by  the  way  similar  to  the  classical  MUSIC  estimator. 
However,  it  is  shown  in  [6]  that  the  accuracy  of  the  linear 
operation  based  methods  [3]-[5]  is  generally  poorer  than  that 
of  the  classical  subspace-based  methods.  On  the  other  hand, 
for  highly  correlated  or  even  coherent  incident  signals,  the 
SUMEW  method  presented  in  [6]  still  relies  on  the  averaging 
techniques.  Although  the  weighted  subspace  fitting  method 
without  eigendecomposition  (WSF-E)  [7]  is  capable  of  re¬ 
solving  the  coherent  signals,  it  still  needs  the  estimate  of  the 


array  covariance  matrix  and  complex  matrix-matrix  products 
to  derive  the  criterion  function,  thereby  indicating  that  the 
WSF-E  method  is  still  computationally  prohibitive 

Recently,  the  methods  termed  ROCK  MUSIC  [8]  and 
ROCKET  algorithm  [9]  based  on  the  MSWF  developed  by 
Goldstein  et  al  [10]  were  proposed  to  high-resolution  spec¬ 
tral  estimation.  Nevertheless,  the  ROCK  MUSIC  technique 
requires  the  forward  and  backward  recursions  of  the  MSWF, 
which  increase  the  complexity  of  the  algorithm.  Moreover, 
the  ROCKET  algorithm  still  needs  complex  matrix-matrix 
products  to  find  the  reduced-rank  data  matrix  and  the  reduced- 
rank  autoregressive  (AR)  weight  vector.  This  implies  that 
additionally  computational  burden  is  included. 

In  this  paper,  we  consider  a  low-complexity  method  of 
weighted  subspace  fitting  (WSF)  for  DOA  estimation.  With 
the  assumption  that  the  training  data  of  one  desired  signal 
are  well  known,  the  novel  method  is  developed.  Firstly,  a 
new  criterion  function  is  derived  based  on  the  MSWF.  The 
novel  signal  subspace  is  obtained  by  the  forward  recursion 
of  the  MSWF,  which  merely  involves  complex  matrix-vector 
products.  And  then,  the  DOA  parameters  can  be  readily 
extracted  by  minimizing  the  novel  criterion  function.  As 
a  results,  the  new  estimator  has  the  attractive  advantages. 
First,  it  does  not  require  the  estimate  of  the  array  covariance 
matrix,  its  eigendecomposition  or  the  backward  recursion  of 
the  MSWF.  Furthermore,  all  operations  included  are  merely 
complex  matrix-vector  products,  thereby  requiring  much  lower 
computational  cost  than  the  existing  subspace  based  methods 
for  direction  estimation.  Second,  all  the  DOA  parameters  of 
the  desired  signal  (with  the  knowledge  of  training  data)  and 
the  ’’interference”  signals  (without  the  knowledge  of  training 
data)  can  be  efficiently  estimated.  Simulation  results  imply  that 
the  proposed  estimator  can  provide  the  comparable  estimation 
accuracy  with  the  classical  WSF  method. 

II.  PROBLEM  FORMULATION 
A.  Data  Model 

Let  us  consider  P  narrow-band  signals  from  distinct 
directions  {0i,  O2,  •  ■  ■ ,  Op}  impinging  upon  a  uniform  linear 
array  (ULA)  composed  of  M  isotropic  elements.  The  M  x  1 
received  noisy  array  data  at  the  fcth  snapshot  can  be  described 
by  the  following  model 

p 

x(fc)  =  ^2  &{0i)si(k)  +  n(fc)  k  =  0, 1,  •  •  • ,  N  —  1(1) 

i- 1 
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where  Si(k)  is  the  scalar  complex  waveform  referred  to  as 
the  ith  signal,  n(fc)  £  CMxl  is  the  additive  noise  vector,  N 
denotes  the  number  of  snapshots,  P  represents  the  number  of 
signals,  and  a(0,)  is  the  steering  vector  of  the  array  toward 
direction  0,  that  is  measured  relative  to  the  normal  of  array, 
and  takes  the  following  form 


a (9i)  = 


Vm  l 


1,  e 


3  Vi 


1  T 


(2) 


where  ipi  =  ‘2lV1  sin  9,  in  which  9i  £  (— n/2,  n/2),  A  and  d 

are  the  wavelength  of  the  carry  signal  and  the  inter-element 
spacing  measured  in  wave  lengths,  respectively,  the  superscript 
(•)T  denotes  the  transpose  operator.  Equation  (1)  can  be 
rewritten  more  compactly  as 


x(fc)  =  A(0)s(k)  +  n(fc)  k  =  0, 1,  •  •  • ,  N  —  1  (3) 


where 


A (9)  =  [a(0!),a (02),  •  •  •  ,a(0P)]  (4) 

s (k)  =  [s1(k),s2(k),---  ,sP{k)]T  (5) 

are  the  M  x  P  steering  matrix  and  the  P  x  1  complex  signal 
vector,  respectively. 

Throughout  this  paper  we  assume  that  M  >  P.  Fur¬ 
thermore,  the  background  noise  uncorrelated  with  signals  is 
modeled  as  a  stationary,  temporally  white,  zero-mean  Gaussian 
random  process,  which  is  also  spatially  white  and  circularly 
symmetric  with  the  second  moments 

E  [n(k)nH(l)~]  =  a\5k^lM  and  E  [n(fc)nT(Z)]  =  0  (6) 

where  Sk,i  is  the  Kronecker  delta  and  I m  denotes  the  M  x  M 
identity  matrix.  We  also  assume  that  all  signals  are  jointly 
stationary,  temporally  white,  zero-mean  complex  Gaussian 
random  processes.  Under  these  assumptions,  the  output  of  the 
array  is  complex  Gaussian  with  zero  mean  and  the  following 
array  covariance  matrix 

Rx  =  E[x(k)xH  (fc)]  =  A(0)RsAh(0)  +  allM  (7) 

where  a2  is  the  noise  variance  and  Rs  =  E[s(k)sH (fc)]  is 
the  signal  covariance  matrix. 


B.  Classical  Weighted  Subspace  Fitting 

Performing  the  eigenvalue  decomposition  of  the  array 
covariance  matrix  Rx  leads  to 


M 

Rx  =  VsAsVf  +  VnAnVp  =  ]T  AiVivf  (8) 

i= 1 

where  Ai  >  >%  >  ■  ■  ■  >  Ap'+i  =  •  •  •  =  A m  =  c£, 
Vs  =  fv-i ,  v2 ,  •  •  ■ ,  vP/]  and  Vn  =  \vP,+1,  vP>+2,  ■  ■  ■ ,  vM],  and 
P'  denotes  the  number  of  uncorrelated  signals.  Accordingly, 
the  eigendecomposition  of  the  sample  covariance  matrix  Rx  = 
jf  J2k=o  x(fc)xff  (k)  is  given  by 


Rx  =  V,A,V 


H 


VnAnV 


H 


M 

J2 Xi ^ 

i= 1 


(9) 


where  Ai  >  A2  >  •  •  •  >  Ap'+i  ~  •  •  •  «  Am  «  0^,  Vs  = 
[v-i,  v2,  ■  ■  • ,  Vp/]  and  Vn  =  |vP/+1 ,  vP<+2 ,  ■  •  • ,  vM]  ■  The  column 
rank  of  Vs  is  in  general  equal  to  the  rank  P'{P'  <  P)  of  the 
signal  covariance  matrix  Rs.  Thus  the  columns  of  Vs  span 
the  P' -dimensional  subspace  of  A(0).  Considering  (7)  and 
(8)  and  performing  some  manipulations,  we  obtain 

Vs  =  A(0)T  (10) 

where  T  €  CPxP  is  the  full  rank  matrix.  Equation  (10) 
forms  a  basis  for  the  classical  signal  subspace  fitting.  9  and  T 
are  unknown  and  can  be  acquired  by  solving  (10).  In  fact,  if 
the  theoretical  Vs  is  replaced  by  the  estimate  Vs,  there  will 
be  no  accurate  solution  to  the  equation  above.  In  this  case, 
one  attempts  to  minimize  some  distance  measure  between  Vs 
and  A(0)T.  For  this  purpose,  the  Frobenius  norm  is  often 
used.  Therefore,  the  SSF  estimator  is  obtained  by  solving  the 
following  non-linear  optimization  problem: 

{0,T}  =  argmin  ||VS  -  A(9)T\\2F.  (11) 

0 .  T 

Since  the  cost  function  above  is  quadratic  with  respect  to 
T,  T  is  easily  obtained.  Inserting  the  least  squares  solution 
T  =  [A(0)ff  A(6*)]_1  Ah (0)Vs  into  (11)  yields  the  following 
equivalent  optimization  problem  without  the  parameter  T: 

Ossf  =  argmm  jtr  (PiVsVf )  j  (12) 

where  =  I m  —  A(9)[A(9)H A(9)]~1  AH (9).  Since  the 
eigenvectors  are  estimated  with  a  quality,  commensurate  with 
the  closeness  of  the  corresponding  eigenvalues  to  the  noise 
variance,  it  is  natural  to  weight  each  eigenvectors  and  lead  to 

9WSf  =  arg  mm  {tr  (piVsWVf )  }  (13) 

where  W  is  the  weighting  matrix  whose  optimal  solution  [2] 
is  given  by  Wopt  =  (As  -  0^1)2^-^ 

C.  Multi-Stage  Wiener  Filter 

The  multi-stage  wiener  filter  presented  by  Goldstein  el  al 
is  to  find  an  approximate  solution  to  the  Wiener-Hopf  equation 
which  does  not  need  the  inverse  of  the  array  covariance  matrix. 
The  MSWF  of  rank  D  based  on  the  data-level  lattice  structure 
is  given  by  the  following  set  of  recursions: 

•  Initialization'.  do(k)  =  si(fc)  and  xo(fc)  =  x(fc). 

•  Forward  Recursion:  For  *  =  1, 2,  •  •  • ,  D: 

hi  =  E[3ci_i(k)d*_1(k)\/\\E[xi_i(k)d*_1(k)}  ||2; 
di(k)  =  hfxi_i(fc); 

Xj(fc)  =  Xi_i(fc)  -  h idi{k). 

•  Backward  Recursion:  For  i  =  D,  D  —  1,  •  •  • ,  1  with 
ep(fc)  =  dD(k): 

Wi  =  E[di^(k)e*(k)\/ E[\ei(k)\2]\ 

£*- i(k)  =  di-i(k)  -  w*£i(k). 

For  the  MSWF  of  rank  M,  the  pre-filtering  matrix  T m  = 
[hi,  h2, . . . ,  hjvr]  is  formed  by  M  matched  filters  in  the 
forward  recursion  of  the  MSWF.  To  reduce  the  computational 
load,  the  MSWF  is  truncated  at  the  Dth  stage  and  the  reduced- 
rank  transformation  matrix  is  Tp>  =  |h| .  h2. ....  h/)].  In  this 
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paper,  we  assume  that  D  >  P'.  Notice  that  the  orthogonal 
matched  filter  h,  £  CM,i  =  1,2,  ...,M  maximizes  the 
real  part  of  the  correlation  between  the  new  desired  signal 
di(k)  =  Xi-\(k)  £  C  at  the  ith  stage  and  the  desired  signal 

di-i(k)  at  the  (i  —  1  )th  stage  subject  to  ||hj(fe)||  =  1,  forcing 
the  desired  signals  between  successive  stages  to  be  in-phase. 
However,  the  blocking  matrix  B,  =1  —  h,h^  guarantees  that 
T m  decorrelates  all  lags  in  the  process  di  ( k )  greater  than  one. 
It  follows  that  the  pre-filtered  covariance  matrix  is  tridiagonal: 


THt)  nr\  _ 


THt) 

D  *Vx0  D 


TH  td  t1 

j~)  -^Vxo  -1-  n' 


Da^x0 

TH  T)  T  rpH  T)  rp 

tXj-iCn  J-  J~)  A.  /lXXo-Ln' 


^xo  D 

5*2 

a3 


J2 

Jd2 


UM 

SM  CTdM  J 


(14) 


=  D 


where  a2d.  =  E  [di(k)d*(k)\,  6i  =  E  [dj(fc)d*_1(A;)] ,  Tn/ 
[hu+ii  hn+i)  ’  ’  ’  i  hM],  and  T M  =  [T x? ,  Tn']- 


III.  LOW-COMPLEXITY  METHOD  FOR  WSF 


A.  Derivation  of  the  Novel  Criterion  Function  for  WSF 

It  is  easy  to  see  from  (14)  that  T|^RXoTd  is  also  a 
tridiagonal  matrix  and  can  be  expressed  as 

TgRX0TD  =  E  [d(fc)dff(fc)]  =  Rd  (15) 

where 


It  readily  follows  from  (14)  that 


TH  T)  rp  _ 

n'  ^x0  J -  D  ~ 


Oix(n-i) 


>D+ 1 


0(M-D)x(D-l)  0(D-l)xl 


(20) 


It  is  shown  in  [10]  that  Xi(fc)  trends  to  become  white  as  the 
stage  of  the  MSWF  increases.  As  a  matter  of  fact,  if  i  £ 
{P',P'  +  1,  •  •  •  ,M  —  1},  Xi(fc)  is  a  temporally  white  random 
process. 

Proposition  1:  Suppose  that  there  are  P'  uncorrelated  sig¬ 
nals  impinging  on  the  ULA,  then  when  the  rank  of  the  MSWF 
is  equal  to  the  number  of  signals,  the  process 

x^  \k)  =  (I  -  hp/hp,)  =  Bp/x^J1(fc)  (21) 


is  a  temporally  white  process,  and  takes  the  following  form 

(22) 

where  the  superscript  (-)^P  1  represents  the  case  where  the 
number  of  incident  signals  equals  P',  the  subscript  (-)p,  refers 
to  the  P'th  stage  of  the  MSWF. 

The  proof  of  Proposition  1  is  seen  in  Appendix  I. 

Corollary  1:  Suppose  that  there  are  P'  uncorrelated  signals 
received  by  the  ULA,  then  when  the  rank  of  the  MSWF  is 
equal  to  or  greater  than  the  number  of  signals,  the  processes 
Xj(fc)  =  B,;Xj_i(fc), i  £  {P',P'  +  -  1}  are 

temporally  white  random  processes. 

Proof:  It  is  easy  to  see  that  the  process  after  the  P'th 
stage  of  the  MSWF  has  the  following  form 


d(fc)  =  Tgx0(fc)  =  [d1(k),d2(k),---,dD(k)]T  (16) 

in  which  dl(k)  is  calculated  by  the  MSWF  algorithm  above. 
For  uncorrelated  signals,  performing  the  eigenvalue  decompo¬ 
sition  of  Rd  leads  to 

D 

TpRXoTp>  =  ES£SE^  +  EnEnE^  =  ^  (17) 

i—l 

where  e^}  ,  i  £  {1, 2,  •  •  • ,  P'}  are  the  largest  eigenpairs  of 
Rd-  Accordingly,  the  eigendecomposition  of  the  estimate  Rd 
is  given  by 

D 

TpRXoTp>  =  ESESE^  +  EnEnE^  =  ^  '  f/jeje^  (18) 

i—l 

where  {?),,  e^}  ,  i  £  {1, 2,  •  •  • ,  P'}  are  the  largest  eigenpairs  of 
Rd.  It  is  easy  to  see  that  the  pre-filtered  matrix  T m  and  its 
estimate  T  M  are  unitary. 

Pre-multiplying  and  post-multiplying  (17)  with  T  /)  and 
Es,  respectively,  and  noting  the  projection  matrix  T  /  , T  = 
I  -  Tn/T^,,  we  obtain 


X,-P  \k) 


B  =  B  iBi-ixjf^fc) 


n  (k) 


n  (k) 


(23) 


where  i  =  P'  + 1,  P'  +  2,  •  •  • ,  M  —  1.  It  follows  from  (22)  and 
(23)  that  the  process  x-P  \k),  i  £  {P',  P'  +  1,  •  •  • ,  M  —  1} 
is  temporally  white.  Thus,  we  obtain  the  Corollary  1.  For 
simplicity,  in  what  follows  we  assume  that  the  number  of 
uncorrelated  signals  is  P'  and  suppress  the  superscript  (■)<l>  1 
of  the  process  x|P  \k),  namely  x,(fc).  ■ 

It  follows  from  Corollary  1  that  when  i  =  {P',  P'  + 
1,  •  •  • ,  M  —  1},  we  have 

Si+i  =  E[di+i(k)d*(k)} 


T/TLEX'E, 

TdT^)RX0Td'E3  (19) 

A(0)RsAH(0)TI?Ei . 


E  [hPixi(fc)xp  i(fc)hi] 
hP  i E  [x^fcjxfL^fc)]  hi 


iH 

Ai+1 


E 


nB' 


A(6>)s  (jfe) 
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+ 


n  BOn(fc) 


\s=i—  1 


H‘ 


(24) 


H 


=  hf+1  £[n(fc)n((f]  []  Bs  h, 

\t=i  /  \s—i—  1  / 

^I-Eh.hf'j  fi-Eh.hf'lh* 


T2hH 

7nrli+l 


=  0. 


Note  that  the  orthogonal  property  of  the  matched  filters  hj  is 
used  in  (24).  So  T  p, RX(|  T u  =  0,  and  (19)  can  be  reduced  to 

T_dEs  (S„  -  all)  =  A(9)RSAH(9)TDES1  (25) 

namely 

TdEs  =  A{9)RSAH(9)TDES  (S8  -  ct2I)_1  .  (26) 

Let 

U  =  TdE  =  Td  [Es  En] 

=  [TdEs  TdE„]  (27) 

=  [Usun] 

where  Us  =  T^ES,  Un  =  T/)En.  Therefore,  (26)  can  be 
reexpressed  as 

Us  =  A(0)K  (28) 

where  K  =  RsAff  (9)TD~ES  (Es  -  oil)'1  e  CPxP' .  It  is 
easy  to  see  that  K  is  of  full  rank.  It  follows  that  Us  spans 
the  signal  subspace.  Thereby,  the  relation  (28)  creates  a  novel 
basis  for  the  SSF,  and  we  have  the  following  new  criterion 
function: 

(0,k)  =  argmin||Us-A(0)K|||.  (29) 

l  J  f?,K 

where  Us  is  the  estimate  of  Us.  Similarly  to  (11),  (29)  is 
also  quadratic  with  respect  to  K.  Thus,  the  parameter  K  can 
be  solved  and  replaced  in  the  criterion  function  above.  For 
the  fixed  unknown  parameter  A (9),  the  solution  for  the  linear 
parameter  K  is 

K  =  At  ( 9 )  Us  (30) 

where  At  ( 9 )  =  [AH  ( 9 )  A  (0)]  1  AH  ( 9 ).  Inserting  (30)  into 
(29),  we  get  the  new  criterion  function  for  SSF  without  K: 

9  =  argmin  ||PiUs|||-  =  argmm  [tr  (piUsUf )  }  . 

(31) 

By  introducing  a  weighting  for  each  signal  vector  estimate  u, , 
i  =  1, 2,  •  •  ■ ,  P',  we  eventually  obtain  the  following  criterion 
function  for  weighted  subspace  fitting: 

9WSf  =  argminllUsW1/2-  A(0)K||f. 

6  ,K 

=  arg  mm  {tr  (piUsWUf )  }  (32) 

where  W  is  the  weighting  matrix  whose  optimal 
value  is  Wopt  =  (Es  —  cr2l)“  in  which 


Es  =  diag{rn,ri2,---  iVp1}- 

It  is  not  difficult  to  show  from  Lemma  1  of  [11]  that  fj, 
and  u,  are  the  Rayleigh-Ritz  fRR)  values  and  RR  vectors  of 
Krylov  subspace  (R.Xo,  rXod0),  which  are  asymptotically 
equivalent  to  the  eigenvalues  and  eigenvectors  of  RXo, 
respectively.  It  follows  that  the  noise  variance  can  be  easily 
estimated  by 


1 


M 


a„  = 


^7  E 


M  -  P 
1 

M  -  1" 
1 

M  -  P> 


i-p'+ i 


p’ 


fr(Rxo)  -  E^ 
1=1 

Pf 

fr(Rxo)  -  E^ 
2=1 


(33) 


However,  the  array  covariance  matrix  RXo  is  unknown 
yet.  Computing  its  estimate  will  increase  the  computational 
cost  of  the  proposed  method.  In  fact,  the  trace  of  RXo 
can  be  computed  by  tr  {rxo{  =  x?x*’  where 

X*  =  XTfj  =  \xi(0) , Xi(l) ,  •  •  • ,  Xi(N  —  1)]T  in  which  fj  = 
[0,  •  •  • ,  0, 1, 0,  •  •  • ,  0]T  and  X  €  CMxN .  It  is  easy  to  see  from 

S- - v- - - ' 

i 

(33)  that  the  array  covariance  matrix  can  avoid  to  be  estimated 
and  the  computational  complexity  of  fr(RXo)  is  reduced  from 
0(M2N)  to  O(MN). 


B.  Computational  Cost  Requirement 

It  should  be  noted  that  the  efficient  implementation  of 
the  MSWF  based  on  the  data-lever  lattice  structure  avoids  the 
formation  of  blocking  matrices,  and  all  the  operations  of  the 
MSWF  only  involve  complex  matrix-vector  products,  thereby 
requiring  the  computational  complexity  of  O(MN)  flops  for 
each  matched  filter  hj,  i  £  {1,  2,  •  •  • ,  D}.  To  fulfil  the  design 
of  the  low-complexity  WSF  estimator,  the  eigendecomposition 
of  the  rank  D  tridiagonal  matrix  is  needed,  which  requires 
0(D3)  complex  product  operations.  Thus,  to  estimate  the 
signal  subspace  matrix  Ts  of  rank  P\  the  computational  cost 
of  the  proposed  method  is  only  0(DM N  +  D2 N  +  D3)  flops. 
However,  the  classical  WSF  estimator  relies  on  the  estimate 
of  the  array  covariance  matrix  and  its  eigendecomposition, 
which  require  0(M2N  +  M3)  flops.  For  the  case  where 
D  «  P'  <C  min(M,  N),  the  computational  complexity  of 
the  proposed  method,  i.e.,  O(DMN)  flops,  is  much  lower 
than  that  of  the  classical  WSF  technique. 


IV.  NUMERICAL  RESULTS 

Example  1  Suppose  there  are  three  uncorrelated  signals 
with  equal  power  in  the  far  field  impinging  upon  a  ULA  with 
16  isotropic  sensors,  whose  spacings  equal  half-wavelength. 
The  true  DOAs  are  { — 5° ,  0° ,  5°} .  The  background  noise  is 
assumed  to  be  a  stationary  Gaussian  white  random  process 
with  zero  mean.  Signal-to-noise  ratio  (SNR)  is  defined  by 
101og((Tg/<T2),  where  a2  is  the  power  of  each  signal  in 
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single  sensor.  The  results  shown  below  are  all  based  on  500 
independent  trials.  The  root-mean-squared  errors  (RMSE’s) 
of  estimated  DOAs  versus  SNR  are  shown  in  Fig.  1.  We 
can  observe  from  Fig.  1  that  the  new  estimator  yields  the 
comparable  estimation  accuracy  with  the  WSF  method  when 
SNR  is  greater  than  -5dB.  The  former  slightly  surpasses  the 
latter  in  estimation  performance  when  SNR  varies  from  -lOdB 
to  -5dB.  AS  SNR  becomes  large,  the  two  estimators  achieve 
the  CRB. 

Fig.  2  shows  the  RMSE’s  of  estimated  DOAs  versus  the 
number  of  snapshots  in  the  case  where  SNR  equals  5dB  and 
the  rank  of  the  MSWF  is  5.  From  Fig.  2,  it  can  be  observed 
that  the  estimation  accuracy  of  the  proposed  method  nearly 
coincides  with  that  of  the  WSF  technique  over  the  range  of 
the  number  of  snapshots  that  we  simulated. 

Example  2  Consider  the  case  where  there  are  three 
signals  impinging  upon  the  ULA  from  the  same  signal  source. 
The  first  is  a  direct-path  signal  and  the  others  refer  to  the 
scaled  and  delayed  replicas  of  the  first  signal  that  represent  the 
multipaths  or  the  ’’smart”  jammers.  The  propagation  constants 
are  {1,  — 0.8  +  j0.6,  — 0.4  +  j0.7}.  The  true  DOAs  are  also 
assumed  to  be  {—5°,  0°,  5°}.  The  results  shown  below  are  all 
based  on  500  independent  trials.  The  RMSE’s  of  estimated 
DOAs  versus  SNR  are  shown  in  Fig.  3.  It  is  easy  to  see  from 
Fig.  3  that  the  estimation  performance  of  the  proposed  method 
is  identical  to  that  of  the  WSF  estimator  when  SNR>-5dB,  and 
the  former  outperforms  the  latter  when  SNR<-5dB.  As  SNR 
becomes  high,  the  RMSE’s  of  the  two  methods  approach  to 
the  CRB. 

For  fixed  SNR  equal  to  5dB  and  the  rank  of  the  MSWF 
equal  to  16,  Fig.  4  indicates  that  the  proposed  method  can 
provide  the  same  estimation  accuracy  as  the  WSF  estimator 
over  the  range  of  the  number  of  snapshots  that  we  simulated. 

V.  CONCLUSION 

An  low-complexity  WSF  method  for  DOA  estimation 
has  been  discussed  in  this  paper.  By  choosing  a  specific 
weighting  matrix,  the  low-complexity  WSF  estimator  yields 
the  comparable  estimation  performance  with  the  classical  WSF 
technique.  Unlike  the  classical  WSF  method,  the  proposed 
method  finds  the  signal  subspace  merely  by  calculating  the 
matched  filters  in  the  forward  recursion  of  the  MSWF,  does 
not  require  the  estimate  of  the  array  covariance  matrix  or  its 
eigendecomposition.  Thus,  the  proposed  estimator  is  compu¬ 
tationally  efficient. 


Therefore,  the  new  observation  data  at  the  first  stage  of  the 
MSWF  is  given  by 


x^(fc)  =  (I  —  hihf )  XQl;(fc) 
a(6»i)aff(6»i) 

l|a(0i)||2 

=  Bm(jfe) 


[a(0i)si(fc)  +  n(fc)] 

(36) 


=  (nB*Jn0fc)- 

Suppose  that  (22)  holds  for  P'  =  K,  namely 

=  (  n  4K\k)  =  (  n  n(fc).  (37) 


Ki=K 


\i—K 


Then,  when  Pf  =  K  +  1,  we  have 


*-(K+ l\k)  = 


U  I  A0 

U=if+1 


x^+1)(fc) 


n  j  x(0  \k)  +  a(0K+i)sK+i{k) 

Ki=K+ 1  / 

/  1  \ 


B 


K+ 1 
1 


IlB 

\i—K 


rW 


i  ^0 


( k ) 


(38) 


IT  B,;  I  a(^tf+i)sif+i(fc)- 


yi=K+i 


Considering  the  orthogonal  property  of  the  matched  filters 
h,,  i  €  {1, 2,  •  •  • ,  K  +  1}  leads  to 


a(0K+i)sK+i{k) 


n  (!-^) 


ti=K+l 
K+ 1 


s(dK+l)sK+l(k) 


i-E1*? 


3{dK+i)sK+i(k)  (39) 


I —  ^K+l^K+i  a( ®K+i  )%:+i( k ) . 


Notice  that  the  orthogonal  matched  filters  h, .  i  € 
{1, 2,  •  •  • ,  A'+l}  are  contained  in  the  signal  subspace,  thereby 
T  k  +  i  T /<- 1  !  is  the  projection  matrix  of  the  column  space  of 
A (6)  =  [a(0i),  a(02),  •  •  ■ ,  &(9k+i)\-  It  follows  that 

(  I]  B  )  a(eK+1)sK+i(k)  =  0.  (40) 

Wif+l  / 


Appendix  I 

The  proof  of  Proposition  1 

When  P'  =  1,  namely  the  case  of  one  signal,  the 
observation  data  reads  as 


Thus,  inserting  (37)  into  (38)  yields 


x 


(tf+l) 

K+l 


(k) 


I]  Bi  n(fc). 


<i=K+ 1 


(41) 


Xo^/c)  =  a(6»i)si(fc)  +  n(k).  (34) 
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